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The recent developments of high intensity uhra-short laser pulses have raised the hopes of observ- 
ing Schwinger vacuum pair production which is one of the important non-perturbative phenomena in 
Quantum electrodynamics (QED). The quantitative analysis of realistic high intensity laser pulses 
is vital for understanding the effect of the field parameters on the momentum spectrum of the 
produced particles. In this study, we analyze chirped laser pulses with a sub-cycle structure, and 
investigate the effects of the chirp parameter on the momentum spectrum of the produced particles. 
The combined effect of the chirp and carrier phase of the laser pulse is also analyzed. These effects 
are qualitatively explained by investigating the turning point structure of the potential within the 
framework of the complex WKB scattering approach to pair production. 
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^ '. INTRODUCTION 

' Schwinger vacuum pair production is one of the important predictions of QED which manifests itself as a non- 
[ perturbative effect of the external field applied to QED vacuum The probability of observing such a phenomenon 

becomes appreciable, when the external field's intensity reaches the scale of the critical field strength E„ = ^ 10^^ 
V/m. However, recent analyses suggest that vacuum pair production may be observed at much lower intensities by 
careful combining and shaping of the laser pulses [3|-|8|. These analyses, together with plans for the experimental 
, realization of field strengths one or two orders below the critical field strength [l^] have raised the hopes of 
^ ' observing such non-perturbative effect ll| . The observation of this effect is crucial for the understanding of other non 



X 



. perturbative phenomena in quantum field theory such as cosmological particle production, 12 1 Unruh and Hawking 



radiation, [ij, [3l heavy ion collisions [3l as well as the Landau-Zener effect [l^ and non-perturbative effects in 
strong field AMO physics [17 



^ ' The quantitative analysis of Schwinger vacuum pair production, for fields varying in time, can be viewed as an 
over-thc-barrier scattering problem in one dimension where the pair production rate corresponds to the reficction 
' coefficient [isl . [Toj . This approach is mathematically equivalent to ca lcu l atin g the particle number expectation in a 



. time dependent background field using the quantum Vlasov equation [20|-|2J| . In a recent study j25j , using the latter 

' approach, the momentum spectrum of the produced particles has been computed for an electric field of the form 

o 

(3 ; E(t) = Eqc ^ cos (tut + if) (1) 



which includes a carrier phase ip in addition to the frequency u and the envelope pulse length t. The momentum 
spectrum was found to be extremely sensitive to these physical pulse parameters, and also to exhibit an interesting 
phase sensitivity between scalar and spinor QED psl [26j . This sensitivity can be understood in terms of quantum 



mechanical resonances, and a quantitative explanation in terms of complex WKB has recently been given in [27|. Here 
we consider a more realistic laser field, incorporating the dependence on the chirp of the laser pulse, which effectively 
includes variations in the frequency, and we investigate the effect of such a chirp on the momentum spectrum. This 
is also important for the study of realistic fields from the experimental point of view; high intensity laser pulses are 
obtained through chirped pulse amplification technique [28j , in which a short laser pulse is chirped and stretched by 
means of a dispersive stretcher, followed by the process of amplification and recompression of the pulse. During this 
process, there might be residual chirp in the field which may introduce frequency variations. Another motivation is 
due to interesting results of chirp effects in strong-field atomic and molecular physics [2^ [s^ . 



Chirped Electric Field 

We consider the electric field E{t) = {0,0, E{t)) as two counter propagating laser beams forming a standing wave 
and the possible effects due to spatial focussing are neglected. Here we analyze the field; 

E{t) ^ Eoer^ cos {bt^ + u}t + (fi) (2) 
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where b is the chirp parameter which causes a hnear variation in the frequency. This type of field configuration is 
represented as E{t) = —A{t) where A{t) is the gauge potential and given as: 

^^^^ = -^^^^^g^e-ra[/(i,^,r,^,6)+g(i,^,r,^,6)] (3) 

where, 

g{t,uj,T,ip,b) = e''^\/j + 26r2(i + 26T2)erfi 

The imaginary error function is defined as: 

erfi(z) = / e-* dt (5) 

When the chirp parameter b is set to zero, these fields reduce to those used in [25|. The effect of the chirp can be clearly 
seen in Fig [T] and [2] The instantaneous frequency of the field linearly increases with the time for a constant chirp 
parameter b. The values of the other parameters such as field intensity Eq and pulse width r arc chosen in the non- 
perturbative regime which may be realized experimentally and given in terms of normalized units (h ~ c ~ m = 1). 
Throughout this paper we use the values of i^o , w and t used in the previous work (25l . [26j , adding various choices 
for the new chirp parameter b. So, we use Eq = O.lEcr, w = 0.05 and r = 100. All the field parameters are given in 
units of the electron mass m. 
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FIG. 1: The thick line shows the electric field ([2]) with chirp parameter & = 0. The dotted [red] line has b — 0.00025, and the 
dashed [gold] curve has b — 0.0005. All plots have Eo = 0.1, cj — 0.05, r = 100 (all the parameters are given in units of the 
electron mass m). Note that the shape of the pulse is very sensitive to the value of the chirp parameter h. 



Computational Formalism 

The calculation of the momentum spectrum of the particles produced by a time dependent electric field can be done 
in several equivalent ways, which we summarize here. The fact that the field is assumed to be spatially uniform means 
that the spatial momentum /c is a good quantum number and we can discuss each mode separately, as the production 
of a particle/anti-particle pair with momenta k and —k. For scalar QED, this involves a mode decomposition of the 
field operator as 



(6) 
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FIG. 2: The thick [blue] hne shows the gauge field Q with chirp parameter b — 0. The dotted [red] line has b = 0.00025, and 
the dashed [gold] curve has b — 0.0005. All plots have Eo — 0.1, cj = 0.05, r — 100 . Note that the shape of the pulse is very 
sensitive to the value of the chirp parameter b. 



Inserting this into the Klein-Gordon equation (— D^I?'' +rn^)^ — 0, yields the following " Schrodinger-like" equation 
for each mode coefficient function u^(t): 

^^+col{t)u^{t)^0 (7) 
where the effective time-dependent frequency of this oscillator problem is 



Wk(t) = y m2 + ^2 ^ _ ^A{t)y (8) 

Here k and k± arc the momenta along and perpendicular to the direction of the electric field, respectivel y. This 
Schrodinger-like problem has scattering boundary conditions, as follows from the usual Bogoliubov expansion 3, 3- 
We make an adiabatic ansatz for ([7]) in terms of the basic leading WKB solutions as: 

u^^t) = ^h^Le-'^^-f*) + -^h^e'Mt^ (9) 

where 



Mt)= / uji,it')dt' . (10) 

We relate the two complex functions ais,{t) and f3k{t) by requiring additionally that Uk(i) follow the constant field 
form: 

Uk(<) = -lio^it) ( ^^e-'^-(*) _ ^^h^e^^A . (11) 
This requires that the Bogoliubov coefficient functions, a\^{i) and /3k (i), satisfy the coupled equations 

«k(0 = l^/^klOe^'^--^*) (12) 

kit) = |^"k(i)e-^'*^« (13) 
Note that these equations preserve the unitarity condition 

|«k(0l'-|/3k(i)|' = l (14) 
which enforces the bosonic field commutation relations. The number of particles produced in the mode k is given by 

AA(k) = |/3k(oo)|2 (15) 
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This particle number ean be expressed in terms of the reflection coefficient for scattering in the Schrodinger-like 
problem ([7]): 



7^(k) = 



/3k(«3) 



Q!k(oo) 



^^'^ - m)--^ (16) 



1+Af{k) ' ' l-7^(k) 



There are many equivalent ways to formulate this scattering computation. In the quantum kinetic equation (QKE) 
approach, one defines a "time-dependent particle number" Af(k, t) = |/3k(0Pj whose time evolution is governed by the 
QKE, inherited from the coupled equations (IT3l) . Alternatively, one can define a time-dependent reflection amplitude 
i?(k, t) = /3k(i)/'^k(i), which satisfies the following Riccati equation [l^, also inherited from the coupled equations 

4caia.(k, t) = (e-^'^^« - i?.^,,,,,(k, Oe^^^-'f*)) , (17) 

where we have added the subscript "scalar" to distinguish from the spinor case below. Using the initial condition 
-Rscaiar(k, -oo) = 0, One evolvcs (H?]) to find 7^scaiar(k) = |i?scaiar(k,oo)p, and hence 7Vscaiar(k) from ([H]). 

For spinor QED, an analogous construction exists, beginning with the Dirac equation instead of the Klein-Gordon 
equation, with the result again formulated either in terms of a QKE or a Riccati equation. Here we use the Riccati 
form, which for spinor QED reads [H,!!!: 



^k(t) ^/rn? + fcj^ {^-2^4,^{t) ^ p2 ^^^2^4>Ut) 
2LJk{t)ik - A{t)) 



^spinor(k, t) — - /,w, _ + -^spinor (k, 0^^*'^'' ) (1^) 



Here (t>-k{t) is the same function as was defined in ([TU|) . Note the changes of sign with respect to the scalar QED Riccati 
equation (jl7p . This reflects the fact that in spinor QED the fields satisfy anti-commutation relations. Computationally, 
the consequence is that the relation between the particle number and the reflection coefficient becomes 



\f . (lA - '^'^P"'°' (^) (iq) 

r "r /vspinor\'^J 



where 72.spinor(k) = |i?spinor(k, oo)p, which is obtained by numerically integrating (|18p . In the subsequent sections, 
we study the particle number spectrum Af(k) for both scalar and spinor QED, obtained by integrating the Riccati 
equations ([T7| and ([T^ . respectively. 

NUMERICAL RESULTS 
Results for zero carrier phase: (fi = 0. 

In this subsection, using the Riccati equation, we compute the momentum spectrum of the produced particles for 
several values of the chirp parameter with the carrier phase ip set to 0. (In all the subsequent computations we take 
k± to be zero.) 

For b ~ 0, the momentum spectrum is centered around origin and oscillations are observed (see Fig. [3]top left). The 
physical origin of the oscillations can be understood as an interference effect between temporally separated reflected 
waves due to multiple bump structure of the scattering potential. The WKB analysis of these effects is explained in 
[27I I in terms of interference between separate complex conjugate pairs of turning points. When the chirp parameter 
b set to 2.5 X 10~*, the momentum spectrum gets shifted in the positive direction and the oscillations become more 
pronounced around bell shaped distribution (see Fig. [3] top right). This result is similar to the effect introduced by 



the carrier phase in [25|. For the higher values of b the spectrum gets shifted in negative direction. The rate gets 
higher compared to the previous cases due to increase in the field's effective frequency. All the spectrums have the 
scalar and spinor results are off-phase due to double valued- ness of the spinor wave function [23|. 

These changes in the momentum spectrum are very sensitive to the value of the chirp parameter. The physical 
explanation for this in the scattering picture is that the corresponding effective scattering potential 

Vk{t) = - {^k - ^A{t)y , (20) 

changes substantially even by small variations of the chirp parameter (see Fig. 2]) 
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FIG. 3: The thick [blue] and dashed [red] curves show the longitudinal momentum spectrum of the produced particles for the 
scalar and spinor cases, respectively. These plots are for zero carrier phase: (p = 0. The other field parameters are Eo — 0.1, 
Lj — 0.05, and r = 100, and from top left to bottom right the values of the chirp parameter are b = 0, b = 0.00025, b = 0.0005, 
and b = 0.00075. Note that, as in the carrier phase fields studied in [2^, the spectra for spinor and scalar QED are out of 
phase. Also notice the strong sensitivity of the spectra to the value of the chirp parameter b. 

Results for non-zero carrier phase: ip — ^. 

To see the effect of carrier phase accompanied with the chirp, we set ip to tt/2 and use the previous values for the 
chirp parameter. When the chirp parameter is set to 0, the scattering barrier becomes completely transparent [i^j for 
various values of momenta which results in almost zero reflection probability (see Fig. [5] top left). When b is set to 
2.5 X 10~^ the momentum spectrum is shifted in negative direction (sec Fig. [S]top ri ght ) . The oscillatory structure 
of the momentum spectrum resembles that of with the parameters 6 = and ip = tt/A [25| . Considering this fact and 
the result represented in Fig. [3] (top right), for relatively small values of &, chirp parameter effectively plays the role 
of carrier phase. However for larger values of b, the oscillations does not have a well defined envelope; there appears 
to be two main peaks in the distribution for scalar case and a single peak for spinor case when b is set to 7.5 x 10~* 
(see Fig. [5] bottom right). 



Turning-point structure 



In this section application of complex WKB to vacuum pair production will be discussed briefly. The results of 
simple gauge fields will be extended by analogy to give a qualitative explanation of the computed momentum spectrum 
of the produced particles pTj . 

As discussed in Section LB, each mode of the field operator satisfies a Schrodingcr-likc ([7]) equation and in the 
scattering picture, the particle number is related to the reflection coefficient, as in (|17p and (fT8)) . Therefore, in this 
approach, complex WKB method Sl'l^- can be used as an analytic way of obtaining the particle number for the 



6 




-20 



FIG. 4: The thick hne shows the scattering potential with 6 = The other values of b are chosen as 0.00025(dotted) and 
0.0005(dashed). The longitudinal momentum is conveniently chosen as — 0.3 to show the effect of chirp parameter on the 
bump structure of the potential. The other parameters are given as £"0 = 0.1, uj = 0.05, r = 100, tp = (upper) and ip — ■k/2 
(lower) 



specified momentum. In obtaining the reflection coefficient, WKB solutions are traced from t — > —00 to < — > 00 to 
find the asymptotic ratio of the WKB coefficients. The dominant contribution to coefficients is obtained when the 
wave function is traced around the points where effective scattering potential Lj{t) is zero. These points are called 
turning points and they appear in complex conjugate pairs when the potential is real. The simplest case assumes that 
the scattering potential just has turning points located on the i mag inary axis in complex conjugate pairs. Previously 
studied gauge potential satisfying this condition is given as 



A{t) = i;ortanh 

where t is the pulse width. WKB result for the particle number density is given as 



K = 



>{t)dt 



(21) 



(22) 



where ti and ^2 are the dominant (ie. closest to the real axis) turning points of the corresponding scattering potential. 
The location of the turning points can be viewed as a function of longitudinal momentum k with r kept fixed, and 
the pair production rate is given as a function of momentum. The momentum spectrum is an exponentially decaying 
distribution and it is in good agreement with the exact result in the quasiclassical regime [2^ . A more general model 
that exhibits the resonance effect has four turning {ti..ti) points where the gauge potential is given as [27[: 



Eo 



Here Eq is the intensity and r represents the width of the field. Pair production rate has been found as [27|, |3; 



N'. 



-2K2 



± 2cos(2a)e' 



-K1-K2 



(23) 



(24) 



where + and — signs represent the scalar and spinor cases respectively. The other quantities are defined as: 
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uj{t) dt 
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uj{t)dt 



<J(K2) 



■Re 



oj(t) dt 



t2 



- ( m ( - 



1 



,1 K 



where the limits of integration correspond to turning points. The result above, involves exponentially decaying terms 
accompanied by an interference term connecting the neighboring turning points (sec Eqn. I24p . The interference term 
is responsible for the oscillations observed in the momentum spectrum. 
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FIG. 5: The thick [blue] and dashed [red] curves show the longitudinal momentum spectrum of the produced particles for the 
scalar and spinor cases, respectively. These plots are for non-zero carrier phase: <^ — n/2. The other field parameters are 
Eo = 0.1, iu — 0.05, and r = 100, and from top left to bottom right the values of the chirp parameter are b — 0, b = 0.00025, 
b = 0.0005, and b = 0.00075. Note that, as in the carrier phase fields studied in [2^, the spectra for spinor and scalar QED are 
out of phase. Also notice the strong sensitivity of the spectra to the value of the chirp parameter b. 

The distribution of the turning points can provide qualitative information about the momentum spectrum. In 
general, the phase integral K is almost linear in t therefore, the dominant contribution to reflection coefficient comes 
from the terms involving turning points which are closest to the real axis. In this model, all the turning points are 
at the same distance to real axis therefore, interference term in (|24p is of the same order of magnitude with the the 
other terms. This fact exp lains why for several values of momenta there can be minima which arc close to zero in the 
momentum spectrum |27| . 

These conclusions can be extended to any field configuration without loss of generality. As for the chirped field, 
there are infinite number of turning points and the relevant ones are shown in Fig. [6] and [7] for various chirp and 
momentum values. For (p ~ and 6 = 0, it is seen that the dominant contribution comes from single central pair 
of turning points. The magnitude of the resonances in the corresponding momentum spectrum (see Fig. [3] top left) 
are small compared to the overall magnitude of the pair production rate. The reason is that interference terms are 
suppressed by the exponential terms and these terms involve turning points which arc further away from the real 
axis than the central pair of turning points (see Fig. [5] top left). Therefore the magnitude of the interference terms 
are small compared to contribution coming from the central pair. On the other hand, for ip = tt/2 and 5 = 0, the 
dominant contribution comes from the two central pairs of turning points (Fig. [7] top left). They are equally far away 
from the real axis as in the model discussed and the argument for four turning point case basically explains why there 
are minima close to zero in the corresponding momentum spectrum (see Fig. [5] top left). 

The effect of the chirp parameter on the spectrum can also be explained by analyzing the distribution of the turning 
points as shown In Fig. [6]and[7l As the chirp parameter changes, the momentum values for which dominating pair of 
turning points gets closest to real axis also change. This results in a shift of the spectrum along the momentum axis. 
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FIG. 6: Contour plots of |cjk(t)P in the complex t plane, showing the location of turning points where ajk(t) = 0. These plots 
are for zero carrier phase: — 0. The other field parameters are Eq — 0.1, lu — 0.05, and r = 100, and from top left to bottom 
right the values of the chirp parameter are b = 0, b — 0.00025, 6 — 0.0005, and 6 = 0.00075, and longitudinal momentum values 
k = 0, k = 0.5, k = —3.0 and k — —3.0, respectively. Notice the strong sensitivity of the locations of the turning points to the 
value of the chirp parameter b. 



The distribution of turning points for non-zero chirp seems to be more complicated; the dominant contribution comes 
from multiple pairs of turning points with multiple interference terms being effective on the momentum spectrum. 
This results in a complicated oscillatory structure. 



CONCLUSIONS 



In conclusion, the two main effects of the chirp parameter on the momentum spectrum can be summarized as 
shifting the spectrum along the momentum axis and increasing the magnitude of the number density at certain 
momenta. The chirp parameter, for relatively small values, has a similar effect on the spectrum as the carrier 
phase. However, when the chirp is increased, the oscillatory behavior of the spectrum gets more complicated. 
The basic features of the computed momentum spectrum such as the shift and the oscillatory structure can 
qualitatively be understood by analyzing the distribution of the turning points. The results by analogy can be 
extended to realistic field configurations. This enables to understand the effect of the important field parame- 
ters such as the carrier phase and the chirp parameter thereby giving useful insight into the structure of external pulse. 

This study was supported by DOE grant DE-FG02-92ER40716. The author thanks G. Dunne for fruitful discussions. 
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FIG. 7: Contour plots of |ajk(t)P in the complex t plane, showing the location of turning points where t<Jk(t) = 0. These plots 
are for non-zero carrier phase: ip — Tr/2. The other field parameters are Eq = 0.1, uj — 0.05, and r — 100, and from top 
left to bottom right the values of the chirp parameter are b = 0, b = 0.00025, b — 0.0005, and b = 0.00075, and longitudinal 
momentum values k = —0.5, k = —3.5, k = —2.0 and k = —0.5, respectively. Notice the strong sensitivity of the locations of 
the turning points to the value of the chirp parameter b. 
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